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ABSTRACT 


A sceady-state  vibrating  object  is  illuminated  with  coherent  light  and 
its  image  is  slightly  misfocused  in  the  film  plane  of  a camera.  The  resulting 
processed  film  is  called  a "time-integrated  specklegram. " When  the  speckle- 
gram  is  Fourier  filtered,  it  exhibits  fringes  depicting  derivatives  of  the 
vibrational  amplitude.  The  direction  of  the  spatial  derivative,  as  well  as 
the  fringe  sensitivity  may  be  easily  and  continuously  varied  during  the 
Fourier  filtering  process.  Besides  the  above  mentioned  advantages,  this 
new  method  is  also  much  less  demanding  than  holographic  interferometry  with 
respect  to  vibration  isolation,  optical  set-up  time,  illuminating  source 
coherence,  required  film  resolution,  etc. 
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Previous  Technical  Reports  to  the  Office  of  Naval  Research 


1.  A.  .r.  Durelli,  "Development  of  Experimental  Stress  Analysis  Methods  to 

Determine  Stresses  and  Strains  in  Solid  Propellant  Grains" — June  1962. 
Developments  in  the  manufacturing  of  grain-propellant  models  are 
reported.  Two  methods  are  given:  a)  cementing  routed  layers  and 

b)  casting. 

2.  A.  J.  Durelli  and  V.  J.  Parks,  "New  Method  to  Determine  Restrained 
Shrinkage  Stresses  in  Propellant  Grain  Models" — October  1962. 

The  birefringence  exhibited  in  the  curing  process  of  a partially 
restrained  p:>olyurethane  rubber  is  used  to  determine  the  stress  associated 
with  restrained  shrinkage  in  models  of  solid  propellant  grains  partially 
bonded  to  the  case. 

3.  A.  J.  Durelli,  "Recent  Advances  in  the  Application  of  Photoelasticity  in 
the  Missile  Industry"--October  1962. 

Two-  and  three-dimensional  photoelastic  analysis  of  grains  loaded  by 
tiressure  and  by  temperature  are  presented.  Some  applications  to  the 
optimization  of  fillet  contours  and  to  the  redesign  of  case  joints  are 
also  included. 

4.  A.  J.  Durelli  and  V.  J.  Parks,  "Experimental  Solution  of  Some  Mixed 
Boundary  Value  Problems" --April  1964. 

Means  of  applying  known  displacements  and  known  stresses  to  the  boundaries 
of  models  used  in  experimental  stress  analysis  are  given.  The  applica- 
tion of  some  of  these  methods  to  the  analysis  of  stresses  in  the  field 
of  solid  propellant  grains  is  illustrated.  The  presence  of  the  "pinching 
effect"  is  discussed. 

6.  A.  J.  Durelli,  "Brief  Review  of  the  State  of  the  Art  and  Expected  Advance 
in  Experimental  Stress  and  Strain  Analysis  of  Solid  Propellant  Grains" — 
April  1964. 

A brief  review  is  made  of  the  state  of  the  experimental  stress  and  strain 
analysis  of  solid  propellant  grains.  A discussion  of  the  prospects  for 
the  next  fifteen  years  is  added. 

6.  A.  J.  Durelli,  "Experimental  Strain  and  Stress  Analysis  of  Solid  Propellant 
Rocket  Motors" --March  1965. 

A review  is  made  of  the  experimental  methods  used  to  strain-analyze  solid 
propellant  rocket  motor  snells  and  grains  when  subjected  to  different 
loading  conditions.  Methods  directed  at  the  determination  of  strains  in 
actual  rockets  are  included. 

7.  L.  Ferrer,  V.  J.  Parks  and  A.  J.  Durelli,  "An  Experimental  Method  to  Analyze 
Gravitational  Stresses  in  Two-Dimensional  Problems"--October  1965. 
Photoelasticity  and  moire  methods  are  used  to  solve  two-dimensional  problems 
in  which  gravity-stresses  are  present. 
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A.  I.  Durelii,  V.  J.  Parks  and  C.  J.  del  Rio,  "Stresses  in  a Square  Slab  ^ 
Bonded  on  One  Face  to  a Rigid  Plate  and  Shrunk"--November  1965.  | 
A square  epoxy  slab  was  Donded  to  a rigid  plate  on  one  of  its  faces  in  ;■ 
the  process  of  curing.  In  the  same  process  the  photoelastic  effects  ' 
associated  with  a state  of  restrained  shrinkage  were  "frozen-in." 

Three-dimensional  photoelasticity  was  used  in  the  analysis. 


9.  A.  J.  Durelii,  V.  J.  Parks  and  C.  J.  del  Rio,  "Experimental  Determination 
of  Stresses  and  Displacements  in  Thick-Wall  Cylinders  of  Complicated 
Shape"--April  1966. 

Photoelasticity  and  moire  are  used  to  analyze  a three-dimensional  rocket 
shape  with  a star  shaped  core  subjected  to  internal  pressure. 


10.  V.  J.  Parks,  A.  J.  Durelii  and  L.  Ferrer,  "Gravitational  Stresses 
Determined  Using  Immersion  Techniques"--July  1966. 

The  methods  presented  in  Technical  Report  No.  7 above  are  extended  to 
three-dimensions.  Immersion  is  used  to  increase  response. 

11.  A.  J.  Durelii  and  V.  J.  Parks.  "Experimental  Stress  Analysis  of  Loaded 
Boundaries  in  Two-Dimensional  Second  Boundary  Value  Problems"-- 
February  1967. 

The  pinching  effect  that  occurs  in  two-dimensional  bonding  problems, 
noted  in  Reports  2 and  4 above,  is  analyzed  in  some  detail. 

12.  A.  J.  Durelii,  V.  J.  Parks,  H.  C.  Feng  and  F.  Chiang,  "Strains  and 
Stresses  in  Matrices  with  Inserts," — May  1967. 

Stresses  and  strains  along  the  interfaces,  and  near  the  fiber  ends,  for 
different  fiber  end  configurations,  are  studied  in  detail. 

13.  A.  J.  Durelii,  V.  J.  Parks  and  S.  Uribe,  "Optimization  of  a Slot  End 
Configuration  in  a Finite  Plate  Subjected  to  Uniformly  Distributed 
Load, "--June  1967. 

Two-dimensional  photoelasticity  was  used  to  study  various  elliptical  ends 
to  a slot,  and  determine  which  would  give  the  lowest  stress  concentration 
for  a load  normal  to  the  slot  length. 

14.  A.  J.  D’orelli,  V.  J.  Parks  and  Han-Chow  Lee,  "Stresses  in  a Split 
Cylinder  Bonded  to  a Case  and  Subjected  to  Restrained  Shrinkage," — 
January  1968. 

A three-dimensional  photoelastic  study  that  describes  a method  and 
shows  results  for  the  stresses  on  the  free  boundaries  and  at  the 
bonded  interface  of  a solid  propellant  rocket . 

15.  A.  J.  Durelii,  "Experimental  Stress  Analysis  Activities  in  Selected 
European  Laboratories"--August  1968. 

This  report  has  been  written  following  a trip  conducted  by  the  author 
tiirough  several  European  countries.  A list  is  given  of  many  of  the 
laboratories  doing  important  experimental  stress  analysis  work  and  of 
the  people  Interested  in  this  kind  of  work.  An  attempt  has  been  made 
to  abstract  the  main  characteristics  of  the  methods  used  in  some  of 
the  countries  visited. 
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16.  V.  J.  Parks,  A.  J.  Durelli  and  L.  Ferrer,  "Constant  Acceleration  Stresses 
in  a Composite  Body"--October  1968. 

Use  of  the  immersion  analogy  to  determine  gravitational  stresses  in 
two-dimensional  bodies  made  of  materials  with  different  properties. 

17.  A.  J.  Durelli,  J.  A.  Clark  and  A.  Kochev,  "Experimental  Analysis  of  High 
Frequency  Stress  Waves  in  a Ring" — October  1968. 

A method  for  the  complete  experimental  determination  of  dynamic  stress 
distributions  in  a ring  is  demonstrated.  Photoelastic  data  is  supple- 
mented by  measurements  with  a capacitance  gage  used  as  a dynamic  lateral 
extensometer . 

18.  J.  A.  Clark  and  A.  J.  Durelli,  "A  Modified  Method  of  Holographic  Inter- 
ferometry for  Static  and  Dynamic  Photoelasticity" — April  1968. 

A simplified  absolute  retardation  approach  to  photoelastic  analysis  is 
described.  Dynamic  isopachics  are  presented. 

19.  J.  A.  Clark  and  A.  J.  Durelli,  "Photoelastic  Analysis  of  Flexural  Waves 
in  a Bar" — May  1969. 

A complete  direct,  full-field  optical  determination  of  dynamic  stress 
distribution  is  illustrated.  The  method  is  applied  to  the  study  of 
flexural  waves  propagating  in  a urethane  rubber  bar.  Results  are 
compared  with  approximate  theories  of  flexural  waves. 

20.  J.  A.  Clark  and  A.  J.  Durelli,  "Optical  Analysis  of  Vibrations  in 
Continuous  Media"--June  1969. 

Optical  methods  of  vibration  analysis  are  described  which  are  independent 
of  assumptions  associated  with  theories  of  wave  propagation.  Methods  are 
illustrated  with  studies  of  transverse  waves  in  prestressed  bars,  snap 
loading  of  bars  and  motion  of  a fluid  surrounding  a vibrating  bar. 

21.  V.  J.  Parks,  A.  J.  Durelli,  K.  Chandrashekhara  and  T.  L.  Chen,  "Stress 
Distribution  Around  a Circular  Bar,  with  Flat  and  Spherical  Ends, 

Embedded  in  a Matrix  in  a Triaxial  Stress  Field" — July  1969. 

A T.hree-dimensional  photoelastic  method  to  determine  stresses  in  composite 
materials  is  applied  to  this  basic  shape.  The  analyses  of  models  with 
different  loads  are  combined  to  obtain  stresses  for  the  triaxial  cases. 

22.  A.  J.  Durelli,  V.  J.  Parks  and  L.  Ferrer,  "Stresses  in  Solid  and  Hollow 
Spheres  Subjected  to  Gravity  or  to  Normal  Surface  Tractions"-- 
October  1969. 

The  method  described  in  Report  No.  10  above  is  applied  to  two  specific 
problems.  An  approach  is  suggested  to  extend  the  solutions  to  a class 
of  surface  traction  problems. 

23.  J.  A.  Clark  and  A.  J.  Durelli,  "Separation  of  Additive  and  Subtractive 
Moire  Patterns"--December  1969. 

A spatial  filtering  technique  for  adding  and  subtracting  images  of  several 
gratings  is  described  and  employed  to  determine  the  whole  field  of 
Cartesian  shears  and  rigid  rotations. 
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24.  R.  J.  Sanford  and  A.  J.  Durelli,  "Interpretation  of  Fringes  in  Stress- 
Holo- Interferometry" --July  1970. 

Errors  associated  with  interpreting  stress-holo-interferometry  patterns 
as  the  superposition  of  isopachics  (with  half  order  fringe  shifts)  and 
isochromatics  are  analyzed  theoretically  and  illustrated  with  computer 
generated  holographic  interference  patterns. 

25.  J.  A.  Clark,  A.  .) . Durelli  and  P.  A.  Laura,  "On  the  Effect  of  Initial 
Stress  on  the  Propagation  of  Flexural  Waves  in  Elastic  Rectangular 
Bars"--December  1970. 

Experimental  analysis  of  the  propagation  of  flexural  waves  in  prismatic, 
elastic  bars  with  and  without  prestressing.  The  effects  of  prestressing 
by  axial  tension,  axial  compression  and  pure  bending  are  illustrated. 

26.  A.  J.  Durelli  and  J.  A.  Clark,  "Experimental  Analysis  of  Stresses  in  a 
Buoy-Cable  System  Using  a Birefringent  Fluid" — February  1971. 

An  extension  of  the  method  of  photoviscous  analysis  is  presented  which 
permits  quantitative  studies  of  strains  associated  with  steady  state 
vibrations  of  immersed  structures.  The  method  is  applied  in  an 
investigation  of  one  form  of  behavior  of  buoy-cable  systems  loaded  by 
the  action  of  surface  waves. 

27.  A.  J.  Durelli  and  T.  L.  Chen,  "Displacements  and  Finite-Strain  Fields  in 
a Sphere  Subjected  to  Large  Deformations" — February  1972. 

Displacements  and  strains  (ranging  from  0.001  to  0.50)  are  determined  in 

a polyurethane  sphere  subjected  to  several  levels  of  diametral  compression. 
A 500  lines-per-inch  grating  was  embedded  in  a meridian  plane  of  the 
sphere  and  moire  effect  produced  with  a non-deformed  master.  The  maximum 
applied  vertical  displacement  reduced  the  diameter  of  the  sphere  by  27 
per  cent. 

28.  A.  J.  Durelli  and  S.  Machida,  "Stresses  and  Strain  in  a Disk  with  Variable 
Modulus  of  Elasticity"--March  1972. 

A transparent  material  with  variable  modulus  of  elasticity  has  been 
manufactured  that  exhibits  good  photoeiastic  properties  and  can  also  be 
strain  analyzed  by  moire.  The  results  obtained  suggests  that  the  stress 
distribution  in  the  homogeneous  disk.  It  also  indicates  that  the  strain 
fields  in  both  cases  are  very  different,  but  that  it  is  possible,  approxi- 
mately, to  obtain  the  stress  field  from  the  strain  field  using  the  value 
of  E at  every  point,  and  Hooke's  law. 

29.  A.  J.  Durelli  and  J.  Buitrago,  "State  of  Stress  and  Strain  in  A Rectangular 
Belt  Pulled  Over  a Cylindrical  Pulley" — June  1972. 

Two-  and  three-dimensional  photoelasticity  as  well  as  electrical  strain 
gages,  dial  gages  and  micrometers  are  used  to  determine  the  stress 
distribution  in  a belt -pulley  system.  Contact  and  tangential  stress  for 
various  contact  angles  and  friction  coefficients  are  given. 
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30.  T.  L.  Chen  and  A.  J.  Durelli,  "Stress  Field  in  a Sphere  Subjected  to 
Large  Deformations"--dune  1972. 

Strain  fields  obtained  in  a sphere  subjected  to  large  diametral  compressions 
from  a previous  paper  were  converted  into  stress  fields  using  two  approaches. 
First,  the  concept  of  strain-energy  function  for  an  isotropic  elastic 
body  was  used.  Then  the  stress  field  was  determined  with  the  Hookean 
type  natural  stress-natural  strain  relation.  The  results  so  obtained 
were  also  compared. 

31.  A.  J.  Durelli,  V.  J.  Parks  and  H.  M.  Hasseem,  "Helices  Under  Load" — 

July  1973. 

Previous  solutions  for  the  case  of  close  coiled  helical  springs  and  for 
helices  made  of  thin  bars  are  extended.  The  complete  solution  is 
presented  in  graphs  for  the  use  of  designers.  The  theoretical  development 
is  correlated  with  experiments. 

32.  T.  L.  Chen  and  A.  J.  Durelli,  "Displacements  and  Finite  Strain  Fields  in 
a Hollow  Sphere  Subjected  to  Large  Elastic  Deformations" — September  1973. 

The  same  methods  described  in  No.  27,  were  applied  to  a hollow  sphere 
with  an  inner  diameter  one  half  the  outer  diameter.  The  hollow  sphere 
was  loaded  up  to  a strain  of  30  per  cent  on  the  meridian  plane  and  a 
reduction  of  the  diameter  by  20  per  cent. 

33.  A.  J.  Durelli,  H.  H.  Hasseem  and  V.  J.  Parks,  "New  Experimental  Method 
in  Three-Dimensional  Elastostatics" — December  1973. 

A new  material  is  reported  which  is  unique  among  three-dimensional 
stress-freezing  materials,  in  that,  in  its  hfated  (or  rubbery)  state 
it  has  a Poisson's  ratio  which  is  appreciably  lower  than  0.5.  For  a 
loaded  model,  made  of  this  material,  the  unique  property  allows  the 
direct  determination  of  stresses  from  strain  measurements  taken  at 
interior  points  in  the  model. 

34.  J.  Wolak  and  V.  J.  Parks,  "Evaluation  of  Large  Strains  in  Industrial 
Applications"--April  1974. 

It  was  shown  that  Mohr's  circle  permits  the  transformation  of  strain  from 
one  axis  of  reference  to  another,  irrespective  of  the  magnitude  of  the 
strain,  and  leads  to  the  evaluation  of  the  principal  strain  components 
from  the  measurement  of  direct  strain  in  three  directions. 

35.  A.  J.  Durelli,  "Experimental  Stress  Analysis  Activities  in  Selected 
European  Laboratories" — April  1975. 

Continuation  of  Report  No.  15  after  a visit  to  Belgium,  Holland,  Germany, 
France,  Turkey,  England  and  Scotland. 

36.  A.  J.  Durelli,  V.  J.  Parks  and  J.  0.  Biihler-Vidal,  "Linear  and  Non-linear 
Elastic  and  Plastic  Strains  in  a Plate  with  a Big  Hole  Loaded  Axially  in 
its  Plane"--July  1975. 

Strain  analysis  of  the  ligament  of  a plate  with  a big  hole  indicates  that 
both  geometric  and  material  non-linearity  may  take  place.  The  strain 
concentration  factor  was  found  to  vary  from  1 to  2 depending  on  the  level 
of  deformation. 
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37.  A.  J.  Durelli,  V.  Paviin,  J.  0.  BUhler-Vidal  and  G.  Ome , "Elastostatics  i 

of  a Cubic  Box  Subjected  to  Concentrated  Loads" — August  1975. 

Analysis  of  experimental  strain,  stress  and  deflection  of  a cubic  box 
subjected  to  concentrated  loads  applied  at  the  center  of  two  opposite 
faces.  The  ratio  between  the  inside  span  and  the  wall  thickness  was 
varied  between  approximately  5 and  121. 

38.  A.  J.  Durelli,  V.  J.  Parks  and  J.  0.  BUhler-Vidal,  "Elastostatics  of 
Cubic  Boxes  Subjected  to  Pressure" — March  1976. 

Experimental  analysis  of  strain^  stress  and  deflections  in  a cubic  box 
subjected  to  either  internal  or  external  pressure.  Inside  span-to-wall 
thickness  ratio  varied  from  5 to  14. 
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Introduc i ion 


Th>’  advantages  common  to  ail  optical  methods  in  experimental  studies 
of  vibrations  ore  that  they  are  noncontacting  and  generally  full-field. 

These  optical  methods  may  be  grouped  into  three  categories  - moir^  methods, 
holographic  interferometric  methods,  and  laser  speckle  techniques. 

A iirief  review  of  the  moird  methods  in  visualization  of  nodes  and 
antinodes  in  vibrating  plate  can  be  found  in  reference  A projected 

ruling  method^^’^^  was  also  developed  for  studying  vibrational  amplitudes. 

A time-dependent  perturbated  ruling  is  observed  when  the  ruling  is  projected 
onto  a vibrating  object.  Photographing  the  perturbated  ruling  with  an  exposure 
time  of  several  vibrational  periods  yields  a time-integrated  moird  fringes 
depicting  tpdal  amplitude  of  the  steady  state  vibration.  These  fringes  can 
also  be  seen  by  eye  in  real-time  if  the  vibrational  frequency  is  approxi- 
mately 20  H or  more.  All  the  moire  methods,  however,  are  limited  to  studies 
of  vibrations  with  relatively  large  amplitudes. 

Powell  and  Stetson  and  other  investigators^ ^ ^ have  demonstrated 
that  time-avera^d  holography  is  well  suitable  for  analysis  of  steady-state 
vibrations  with  very  small  amplitudes.  In  the  reconstruction  of  the  holo- 
gram which  recorded  a steady-state  vibrating  object,  a time-averaged  fringe 
pattern  is  observed  which  is  described  by  a zero-order  Bessell  function 
containing  the  vibrational  amplitude  in  its  argument.  Nodal  areas  are  also 
clearly  indicated  by  the  regions  of  the  brightest  reconstruction.  An 
extension  of  this  technique,  the  modulated  holography ^ ^ \ can  be  used  to 


reduce  the  sensitive  i; 


•jnd  to  determine  the  relative  phase.  Shortcomings  of 


these  tediniques,  however,  are  Chose  associated  with  holographic  interfero- 
metry, such  as  stringent  requirement  on  mechanical  and  thermal  stabilities, 
Che  fringe  localization  problem,  etc. 

The  speckle  pattern  produced  by  the  random  scattering  of  a diffused 
object  when  illuminated  by  coherent  light  acts  as  an  information  carrier 
about  the  object  surface.  By  simply  watching  the  speckle  pattern  on  a 
surface  wtiile  it  is  undergoing  ste  .dy  st.i.ce  oscillation,  one  can  observe 
nodal  patterns  in  real  time^^^’^^’^  This  provides  a quick  and  simple  way 
of  scanning  for  the  appropriate  frequencies  to  be  studied  in  greater  detail. 
Butters  and  Leendertz  demonstrated  that  a time-averaged  speckled  fringe 
pattern  depicting  vibrational  amplitudes  similar  to  that  of  time  averaged 
holography  could  be  obtained.  The  sensitivity  of  their  technique  i':  abo-t 


the  same  order  as  that  of  holography.  Nevertheless,  it  is  less  demanding 
on  film  resolution  and  speckle  fringes  localize  on  the  object  surface. 


Anothr-r  virtue  of  speckle  interferometry  over  holographic  interferometry  is 

(14,15,16) 

tha^  it  permits  in-plane  displacements  to  be  measured  independently 

...  , . (17,18) 

This  advaritage  is  particularly  useful  for  in-plane  vibration  analysis 


Another  interesting  technique  for  vibration  studies  is  the  speckle- 

( 19 ) 

shearing  interferometer  introduced  by  Hung  and  Tavior  . The  technique 

measures  derivatives  of  vibrational  amplitudes  directly;  thus  it  greatly 
facilitates  the  analysis  of  stresses  resulting  from  vibrations.  An  improved 
•/"rsion  was  reported  by  Hung,  Rowlands  and  Daniel^  \ In  their  method,  the 
object  under  study  was  illuminated  by  a beam  of  coherent  light  and  was 
imag«*d  by  a camera  iiaving  a mask  with  four  apertures  in  front  of  the  lens. 

A photographic  plate  was  exposed  in  the  misfocused  image  plane  while  the 
object  was  undergoing  steady-state  vibration.  The  purpose  of  misfocusing 


2 
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was  to  produce  a sh'-arinf'.  effect.  The  proces.sed  plate  was  Fourier  filtered 
to  yield  four  families  of  fringes,  each  depicting  the  derivatives  of  the 
vibrational  amplitudes  with  respect  to  the  direction  parallel  to  the  line 
joining  the  centers  of  each  aperture  paii  . The  sensitivity  of  the  method 
was  r,.ainly  controlled  bv  the  amount  of  shearing.  However,  both  the  deri- 
vative directions  as  well  as  the  sensitivity  were  fixed  after  recording. 

The  present  paper  describes  a new  optical  technique  whereby  the  deri- 
vatives of  vibrational  amplitudes  with  respect  to  all  directions  can  be 
recorded  at  one  time.  Moreover,  the  sensitivity  can  be  varied  after 
recording . 

Tne  fact  the  technique  measures  derivatives  of  vibrational  amplitudes 
is  particularly  useful  for  studying  flexural  stresses  induced  by  vibrations. 
Since  strains  due  to  vibrational  bending  are  proportional  to  the  second  deri- 
vatives of  the  amplitudes,  it  therefore  only  requires  one  differentiation 
to  obtain  strains  as  in  contrast  to  those  amplitude  measuring  techniques 
which  need  to  differentiate  the  data  twice. 

The  examples  presented  in  Fig.  4 were  prepared  at  low  levels  of 
sensitivity  in  order  to  show  the  principle  of  this  technique.  For  actual 
quantitative  studies,  the  systems  should  be  adjusted  to  yield  higher  fringe 
densities  in  order  to  facilitate  accurate  determination  of  the  second 
derivatives . 


Description  of  the  Method 

Figure  1 illustrates  the  schematic  diagram  for  recording.  The  object 
under  study  is  illuminated  by  a collimated  beam  of  coherent  light  inclined  at 
an  angle  0 to  the  viewing  direction.  It  is  imaged  by  a camera  with  the 
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aperture  fully  opened,  A phototnraphic  plate  is  placed  at  a small  distance  S 
from  the  image  plane.  While  the  object  is  undergoing  steady  state  oscillations 
the  photographic  plate  is  exposed  with  an  exposure  time  equal  to  many 
vibrational  periods.  The  processed  photographic  plate  is  referred  to  as  a 
"time  integrated  specklegram . " It  will  be  shown  in  the  following  that  a 
i rin.’e  pat'ern  con-  ur  ing  ’ h<;  derivatives  of  the  vibrational  amplitudes  with 
r.,v;;.,-t  I'  any  a-sirod  ' tian  can  be  extracted  from  the  specklegram  with 
vari  idl  .n-nsitiv  I f,-  '’V  ^ .■ourl'-r  filtering  technique. 


Ana  Lyr,  is  of  Recording  ! rocess 


The  full  aperture  of  the  lens  can  be  considered  as  an  infinitely  dense 
collection  of  small  aperture  pairs.  Let  us  focus  on  the  function  of  one 
aperture  pair.  Let  the  line  joining  the  centers  of  the  aperture  pair  be 
parallel  to  the  r-direction  as  shown  in  Figure  2.  and  d be  their  separation. 
The  two  images  of  the  structure,  one  focused  by  each  aperture,  are  coincident 
in  the  image  plane.  Due  to  the  film  plane  being  misfocused,  these  two 
images  are  sheared  laterally  in  the  r-direction.  The  amount  of  the  shearing 
^r'  in  the  film  plane  is 


6r ' = 


(1) 


where  D.  is  the  distance  of  the  image  plane  from  the  lens,  and  S is  the 
1 

misfocus.  The  equivalent  shearing  6r  in  the  object  plane  is 


6r 


(2) 


where  m is  the  magnif ication . 

Obviously  this  arrangement  brings  the  ray  scattered  from  a point  P(x,y)  to 
meet  the  ray  scattered  from  a neighboring  point  P(x+6x,  y+6y)  in  the  film 
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plane  where  5x  = 6r  coi:,),  and  6y  = 6r  sin(};.  If  the  two  aperfores  are  small, 
they  can  be  considered  as  a two-point  source,  and  the  resulting  interference 

is  a ruling  whose  lines  are  perpendicular  to  the  r-direction.  The  frequency 

, . . . , (15,1b) 

f'  of  the  ruling  is  given  by 


A(D^+S) 


where  X = wavelength  of  the  laser  light. 

Hence  the  equivalent  frequency  f in  the  object  is: 


xToTTsT 


m^D. 

1 

\s  (D^+S) 


The  ruling  is  distorted  by  the  randomness  of  the  speckle  pattern.  It 

can  be  shown  that  the  intensity  distribution  I(x,y)  of  the  speckle  ruling  is 

^ , (19,20) 

represented  by 

I(x,y)  = a^[l  + cos(27rfr  + x)3*  (5) 

where  x ~ random  speckle  phase 
a - light  amplitude 
X - rcos4i;y  = rsin()) 

It  has  been  assumed  in  the  above  analysis  that  the  object  is  in  its  stationary 
position.  Suppose  that  the  object  is  undergoing  steady  state  sinusoidal 
vibration  with  a displacement  only  in  the  z-direction  represented  by: 

W(x,y,t)  = A(x,y)  cos(ut+a)  (6) 


’’Por  simplicity,  the  magnification  is  assumed  to  be  unity  and  so  the  image  plane 
would  take  the  same  coordinates  as  the  object  plane  except  it  is  inverted. 

This  assumption  is  carried  on  hereafter. 
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where 


W(x,y,t)  is  the  displacement  in  the  z-direction 
A(x,y)  - amplitude;  function 


u)  - circular  frequency 

a - arbitrary  phase 

The  relative  displacement  6W(t)  between  the  point  P(x,y)  and  the  neighboring 
point  P(x+6x,  y+6y)  is  thus: 

6W(x,v,t)  = W(x+6x,yt6y ,t)  - W(x,y,t)  (7) 

This  relative  displacement  will  cause  a relative  optical  path  change  and  hence 

the  phase  between  the  light  scattered  from  the  two  points.  With  the  assumption 

that  the  object  size  is  small  compared  to  D^,  this  relative  phase  A change  can 

be  shown  to  be  approximately  given  by: 

A = (1+COS0)  • 6W(x,y,t)  fS) 

4 

which  results  in  a local  optical  shift  of  the  speckle-ruling.  Thus  equation 
(5)  in  which  I(x,y)  is  now  time  dependent  becomes; 

I(x,y,t)  = a^[l  + cos(  2irfr+x+A)]  (9) 

The  exposure  E(x,y)  on  the  film  for  a number  of  vibrational  period  is  there- 
fore an  integration  of  the  intensity  over  the  exposure  time  T given  by: 

E(x,y)  = k I(x,y,t)dt  (10) 

o 

where  k is  the  photographic  film  constant 
The  al-)Ove  integration  yields: 

E(x,y)  = kTa^fl  + cos(2iTfr+x)  ' (H) 

where 

6 = ^^1  + cos0)[A(x+6x,y+6y ) - A(x,y)]  (12) 

Equation  (11)  represents  a carrier  of  approximately  constant  frequency 
modulated  by  a Bessel  function  of  zero-order.  Nulling  of  the  carrier 
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r 1 

i 
1 
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which  will  be  Identifi''d  a.-:  moire  fringes  occurs  when  0 equals  to  the  null  ^ 

i 

roots  of  the  zero-order  Bessel  function.  Equation  (12)  may  be  rewritten  as:  1 

2Tt,.  . r A(x+6x,yt6y)  - A(x,y)-, 

6 = ;p(i  + COS0)  6r[ (3^3) 

2ti,,  r ri5A(x,y)-, 

= + cosO)  6r[ ^ — J 

X or 

If  6r,  the  lateral  shear  is  small,  (|^)  approaches  (|^)  • Hence  the 
fringe  pattern  depicted  by  equation  (11)  is  a contour  map  of  the  derivatives 
of  the  vibrational  amplitudes  with  respect  to  the  r-direction. 

It  is  seen  in  equation  (13)  that  6r,  the  amount  of  shearing,  is  a 
sensitivity  factor,  and  it  depends  upon  d,  the  separation  of  the  aperture 
centers.  Also  6r  is  related  to  the  frequency  of  ruling  by  equation  (4).  Thus 
an  aperture- pair  with  greater  separation  produces  a greater  shearing  and  hence 
greater  carrie’"  frequency.  The  direction  of  the  shearing  is  parallel  to  the 
line  -joining  the  aperture-pair  centers  i.e.,  it  is  perpendicular  to  the  lines 
of  the  ruling.  With  the  orientation  of  the  aperture-pair  varied,  both  the 
shearing  direction  and  the  ruling  orientation  are  changed  accordingly. 

We  have  just  shown  that  an  aperture-pair  is  indeed  an  interferometer 
measuring  the  derivatives  of  the  vibrational  amplitudes  with  respect  to  a 
direction  parallel  to  the  line  joining  the  aperture-pair  centers  and  with 
sensitivity  depending  upon  the  aperture  separation.  This  information  is 
represented  by  a family  of  fringes  formed  by  moire  phenomenon  of  two 
slightly  mismatched  speckle-rulings  interpreted  by  equation  (11).  Since 
the  whole  lens  aperture  consists  of  many  aperture  pairs  with  aperture 
separation  varying  from  0 to  R,  the  diameter  of  the  lens,  and  with  the 
orientation  varying  from  (fi  = 0 to  <fi  = 360°,  recorded  on  the  specklegram 
are  rulings  of  all  orientations  and  of  frequency  varying  from  0 to 
R/X(D.+S).  Each  ruling  forms  a fringe  family  of  its  own.  Therefore  a 
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specklef^ram  containu  a :.uperposition  of  a multiplicity  of  moire  fringe 
families.  Eacli  fringe  family  depicts  vibrational  amplitudes  with  respect 
to  a direction  perpendicular  to  the  ruling-line  and  with  sensitivity  depends 
on  the  ruling  frequency.  Thus  the  derivatives  of  the  vibrational  amplitudes 
with  respect  to  all  directions  are  recorded  in  the  specklegram  with  a 
variable  range  of  sensitivity  after  recording  the  upper  limit  depending  on 
the  highest  frequency  which  the  recording  lens  is  capable  of  passing. 

Readout  Process  by  t^ourier  Filtering 

To  readout  the  fringe  family  depicting  amplitude-derivative  with  respect 
to  a particular  direction  and  with  a desired  sensitivity,  a readout  device 
shown  in  figure  3 is  employed.  A beam  from  a laser  is  expanded  by  a micro- 
scopic objective  and  is  spatially  filtered  by  the  pinhole.  The  beam  diverged 
from  the  pin  hole  is  focused  by  a lens  to  a point  lying  in  a plane  known  as 
the  Fourier  filtering  plane.  The  time-integrated  specklegram  is  placed  in 
the  Input  plane  as  shown  and  is  imaged  by  a camera  behind  the  Fourier 
filtering  plane.  It  will  be  shown  that  by  inserting  a small  aperture  off-center 
in  the  q-direction,  a fringe  pattern  depicting  the  vibrational  amplitude 
derivative  with  respect  to  the  n-direction  is  displayed  in  the  image  plane  of 
the  camera  (output  plane)  with  a sensitivity  proportional  to  the  distance 
of  the  aperture  from  the  optical  center. 

It  can  be  shown  that  what  appears  in  the  Fourier  filtering  plane  of 

Figure  3 is  the  Fourier  transform  of  the  amplitude  transmittance  of  the 

(21) 

specklegram  multiplied  by  a quadratic  phase  term.  This  represents  a 

transformation  from  a spatial  domain  to  a frequency  domain.  In  other  words, 
rulings  of  different  frequency  are  separated  in  the  Fourier  filtering  plane. 

If  a small  aperture  is  located  at  a distance  q from  the  optical  axis  in 
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the  n-diroction,  a rulln;’  orthogonal  to  the  n-direction  and  of  frequency 
f is  given  by 

f = V (13) 

X 'Z 

where  X = wavelength  of  light  used  for  filtering  is  passed.  As  analyzed 
above  this  ruling  corresponds  to  a shearing  6r  in  the  p-direction  given  by: 


6r 


XS(D.+S) 

= a i 

X ' Z D . 


(14) 


However,  areas  corresponding  to  the  nulling  of  the  carrier  as  described  in 
equation  (11)  makes  no  contribution.  Thus  those  areas  will  appear  dark  in 
the  output  plane.  Therefore  a fringe  pattern  is  observed  on  the  image  that 
depicts  the  amplitude-derivative  with  respect  to  the  n-direction  accompanied 
by  a sensitivity  factor  of  6r  given  in  equation  (l5)-  By  varying  the 
n-direction,  any  desired  derivative  direction  can  be  chosen  with  an  adjustable 
sensitivity  depending  on  the  distance  q. 


Experiments 

A rectangular  plate  of  dimensions  3.75  x 5.0  cm  clamped  along  its  boundaries 
was  chosen  for  demonstration.  The  plate  was  driven  by  an  audio  speaker  so  that 
it  was  vibrating  at  its  fundamental  frequency.  The  set  up  shown  in  Figure  1 
was  used  for  recording  whereby  a spectra  physics  125  laser  was  used  for  the 
illuminating  source.  The  processed  time-integrated  specklegram  was  Fourier 
filtered  with  the  scheme  of  Figure  3.  The  specklegrams  were  recorded  on  Agfa 
10E75  film,  typically  using  exposure  times  of  1-5  seconds.  Photographic 
processing  included  5 minutes  in  Kodak  D-19 , 1 minute  stop,  5 minutes  Fix  and 
10  minutes  wash.  Figure  4 shows  the  fringe  patterns  extracted  from  the 
specklegram  depicting  the  derivatives  of  the  vibrational  amplitudes  with 
respect  to  the  various  directions  as  shown,  and  with  different  sensitivities. 
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Summary  and  Conclusion 


i 

t 

[ 

i_ 

i A optical  technique}  has  been  described  which  permits  derivatives  of 

vibrational  atnplitades  to  be  measured  directly.  One  outstanding  feature  of 
the  technique  is  that  it  provides  variability  after  recording,  both  in 
derivative  direction  as  well  as  sensitivity.  The  method  can  be  viewed  as 
a holographic  interferometric  method  in  which  every  point  on  the  object 
acts  as  a reference  for  all  its  neighboring  points.  However,  several 
limitations  associated  with  holography  have  been  alleviated.  These  include; 

(1)  simple  experimental  set  up;  (2)  less  demanding  on  stabilities,  (3)  relax- 
ation of  film  resolution  requirement  and  (4)  providing  greater  range  of 
sensitivity.  A shortcoming  of  the  technique  is  that  fringe  visibility  is 
generally  poor. 

It  should  be  noted  that  the  mechanical  shift  of  speckle-rulings  due  to 
lateral  displacement  (displacement  orthogonal  to  line  of  sight)  is  assumed 
negligible.  Should  the  lateral  displacement  be  appreciable,  the  mechanical 
shift  must  be  considered.  The  combined  effect  of  the  optical  and  mechanical 
shift  of  speckle-rulings  will  be  analyzed  in  a future  paper. 
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FIGURE  : FUNCTION  OF  AN  APERTURE- PAIR 


FILTERING  APERTURE 


FIGURE  3:  FOURIER  FICTIRING  SCHEME  FOR  EXTRACTING 

DISPLACEMENT-DERIVAHVE  FRINGES 
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A Steady-State  vibrating  object  is  illuminated  with  coherent  light  and 
its  image  is  slightly  misfocused  in  the  film  plane  of  a camera.  The  resulting} 
processed  film  is  called  a "time- integrated  specklegram.  When  the  speckle- 
gram  is  Fourier  filtered,  it  exhibits  fringes  depicting  derivatives  of  the 
vibrational  amplitude.  The  direction  of  the  spatial  derivative,  as  well  as 
the  fringe  sensitivity  may  be  easily  and  continuously  varied  during  the 
Fourier  filtering  process.  Besides  the  above  mentioned  advantages,  this 
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